Incommensurate phases in ferromagnetic spin-chains with weak antiferromagnetic 
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We study planar ferromagnetic spin-chain systems with weak antiferromagnetic inter-chain interac- 
tion and dipole-dipole interaction. The ground state depends sensitively on the relative strengths of 
antiferromagnetic exchange and dipole energies k = J'a 2 c/{g L HB) 2 - For increasing values of k, the 
ground state changes from a ferromagnetic via a collinear antiferromagnetic and an incommensu- 
rate phase to a 120° structure for very large antiferromagnetic energy. Investigation of the magnetic 
phase diagram of the collinear phase, as realized in CsNiF3, shows that the structure of the spin 
order depends sensitivly on the direction of the magnetic field in the hexagonal plane. For certain 
angular domains of the field incommensurate phases appear which are separated by commensurate 
phases. When rotating the field, the wave vector characterizing the structure changes continuously 
in the incommensurate phase, whereas in the commensurate phase the wave vector is locked to a 
fixed value describing a two-sublattice structure. This is a result of the competition between the 
exchange and the dipole-dipole interaction. 



I. INTRODUCTION 

Ferromagnetic spin chains have been studied exten- 
sively, experimentally as well as theoreticallytrtj. Typical 
systems are CsNiF3 and RbFeCl3 which show quasi one- 
dimensional behavior due to the small lattice constant 
along the c-axis. Nearest neighbor spins couple with a 
ferromagnetic exchange interaction along the spin chain. 
A planar anisotropy of the same order of magnitude is 
found because of the non-vanishing orbital angular mo- 
mentum of the magnetic ions. Perpendicular to the chain 
axis the magnetic ions are located on a triangular (hexag- 
onal) lattice with a much larger lattice constant a, i.e. the 
magnetic ions form a simple hexagonal lattice. Thus, at 
high temperatures a pronounced one-dimensional behav- 
ior can be found experimentally. Planar ferromagnets 
still have rotational symmetry, so that no long range 
order exisjte. Nevertheless, spin waves were measured 
for q ^ 0q At low temperatures the systems undergo 
a phase transition to a three-dimensional ordered phase 
for which the antiferromagnetic exchange interaction and 
the dipole-dipole interaction are responsible. The precise 
structure of the order in the hexagonal plane depends 
sensitively on the relative strengths of these two compet- 
ing interactions. RbFeCi3 has a nearly 120° structure be- 
cause the exchange energy is large. In contrast, CsNiF 3 
has a collinear orientation of the spins due to the large 
dipole energy. Up to now the three-dimensional ordering 
has not attracted much interest, but interesting behavior 
is expected due to the competing interactions which can 
lead to frustration. _ . 

Recently Baehr et al£3 have measured the magnetic 
excitations in the three-dimensional ordered CSN1F3 
(T < TV = 2.7 K) by inelastic neutron scattering. For 



the first time the magnitude of the antiferromagnetic ex- 
change energy in the plane has been determined by mea- 
suring the dispersion relation. The important role of the 
dipole-dipole interaction is manifested by the disconti- 
nuity in the magnon energy at the zone center (q=0). 
This is an effect of the long range interaction. In this 
paper we want to study the whole range from vanish- 
ing exchange energy up to the region where the dipole 
energy is small. It will be shown that instead of the frus- 
trated 120° structure an incommensurate phase will be 
established. The dispersion relations are evaluated for 
the commensurate phases. In a second step we investi- 
gate the magnetic phase diagram for the collinear phase, 
especially for CSN1F3, and calculate the critical fields. 
The Ni atoms form a simple hexagonal lattice with lat- 
tice constants a jy b = 6.2lA and c = 5.2A the Lande 
factor gL = 2.250 and the spin S = 1. The outline of 
the paper is as follows. In chapter II we introduce the 
model, in chapter III the ground state for vanishing field 
strength is determined and in chapter IV we study the 
excitations of the ferromagnetic and the collinear phase. 
In chapter V the collinear antiferromagnetic phase with 
finite fields in the hexagonal plane is investigated. In the 
appendix we summarize the main relations for the dipole 
energy. 



II. MODEL 

The starting point of our investigation is the Heisen- 
berg Hamiltonian 



1 



a,/3 1,1' 



/ 

I 



(1) 



Here J denotes the ferromagnetic nearest-neighbor in- 
trachain interaction, A the single-ion anisotropy, J' w the 
the dipole-dipole interaction 
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i indicates positions on one and the same spin chain, 
whereas I indicates all spin positions. The notation for 
xi, used here, is given in the appendix. Ho is an external 
field perpendicular to the chain axis. We are interested 
in the case of a planar ferromagnetic chain, i.e. J, A > 0, 
which means that the spins are forced to lie in the hexag- 
onal plane perpendicular to the chain. The following cal- 
culations are valid only for weak inter-chain interactions, 
i.e. J, A ^> J', (gLfiB) 2 /a 2 c, where a indicates the lattice 
constant of the triangular lattice in the plane and c the 
lattice constant of the chain. Because of the long-range 
of the dipolar interaction, summations over the whole 
system arc performed best by means of the Ewald sum- 
mation techniquelld (s. appendix). In general, we can 
differentiate between one-dimensional behaviour at high 
temperatures due to the large intra-chain energy and the 
three-dimensional properties at low temperatures. 
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The classical ground state is determined by the maximal 
eigenvalue of 
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with the wave- vector q constrained to the plane (q z = 0). 
The eigenvalues of the three by three matrix (Eq. ^) are 
given by 

Ai = V*' 



= ^ [V™ + V™ ± V(KT - W + 4 « V f ) - (8) 

because the off-diagonal components A^ z and A^ z vanish 
for wave vectors in the plane. An out-of-plane orientation 
of the spins increases the energy (V^ z ) because of the pla- 
nar anisotropy. Thus, the wave vector, which maximizes 
A2, describes the classical ground state. In Fig. ([j]) the 
exchange energy J q and the the three components of the 
symmetrical dipole tensor A^ 3 are plotted (a/c = 2.39 
and J' = 1). The maximal value for pure antiferromag- 
netic exchange is found to be for qo = 47r/3a(l, 0, 0), i.e. 
a 120° structure. 



III. GROUND STATES FOR H = 

First we study the ground state of the three- 
dimensional system without an external field. Due to 
the pronounced one-dimensional behaviour the classi- 
cal ground state consists of ferromagnetic spin-chains 
with the spins ordered in the plane perpendicular to the 
chains. The spin structure within the hexagonal plane 
then depends only on the antiferromagnetic exchange and 
the dipolar interaction. For the following investigation 
we introduce a dimensionless parameter k, which is the 
ratio of the antiferromagnetic exchange to the dipolar 
interaction (KcsNiF 3 = 0.79): 
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Fourier transformation of the Hamiltonian (Eq. (Q)) 
yields 

H = - £ £ ( J q <5 Q/3 - A8 az 8 Pz + J q (5 Q/3 + Af) S^Si^ , 

a, (3 q 

(4) 

with the nearest-neighbor exchange energies (intrachain 
and interchain, J, J' > 0) 




FIG. 1. Antiferromagnetic exchange energy Jq and dipolar 
interaction Aq 13 for selected directions in the Brillouin zone. 
The solid line represents the isotropic exchange, the dashed 



line Aq y , the dot-dashed line A^ x 



and the dotted line A^ y . 



The dipole components are given in units (<?lMb) / a 3 an d 
the exchange in units J'. 

For pure dipolar interaction the minimal value is at 
wave vector q = 0, a ferromagnetic orientation in corre- 
spondance to the two-dimensional counterpart, the tri- 
angular latticeO. The ferromagnetic structure is stable 
as long as the antiferromagnetic exchange energy is lower 
than the dipolar interaction leading to 
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2 7 r/v / 3a(0,l,0). (9) 



Due to the semiconvergence of the dipole sums for a fer- 
romagnetic structure (q = 0), the dipole components are 
not analytic at the zone center. For spherically shaped 
systems the values for the dipole tensor are given in the 
appendix. For CSN1F3 Eq. (m reduces to the inequality 



J' < 7 mK 



K < K c l = 0.20 



(10) 



As a result of the hexagonal symmetry the ferromagnetic 
ground state is continuously degenerate. 

When the exchange energy exceeds the above inequal- 
ity (Eq. (0)) the classical ground state chajiges to a 
collinear antiferromagnetic spin orientationli3. There 
the spins are orientated along the lattice axes so that 
there exist only three discrete states called domains 
(Fig. ||). Domain A can be described by the wave vec- 
tor qi, domains C and B by q2 = 7r(l, 1/v3, 0) and 
q 3 = 7r(l, — 1/\/3, 0) respectively. These wave vectors 
correspond to points P, X and X' in the Brillouin zone 
(Fig- 1). 




FIG. 2. The Brillouin zones of the hexagonal plane. The 
hexagon is the crystallographic and the rectangle (dashed) is 
the magnetic one. Included are the three directions for the 
scattering measurement which are realized in the domains. 




FIG. 3. Transition from the collinear antiferromagnetic 
phase, points P, X and X' , towards the three-sublattice phase 
J. The intermediate phases are incommensurate, i.e. the 
wave vectors, defining the structure, have no commensurate 
relation to the periodicity of the underlying lattice. 

The incommensurate wave vector behaves near P as 
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i.e. a square root dependence. For increasing exchange 
energy point J is reached asymptotically (s. Fig. 




For increasing exchange energy the collinear ground 
state does not change to a frustrated 120° structure as 
one might have assumed. Rather, the system changes to 
an incommensurate phase. This results from the fact that 
the dipole energy depends linearly on the wave vector 
near point J, in contrast to a quadratic behavior for the 
exchange energy (Fig. Q). Therefore, the wave vector 
for the lowest ground state energy changes continuously 
from point J towards point J, which would correspond 
to a three sublattice state. Point J is reached only for 
vanishing dipole energy. Due to the hexagonal symmetry 
the same situation holds for points X and X' ', i.e. the 
three-domain structure from the collinear phase survives 
in the incommensurate regionHQ. The three domains and 
the route to the three sublattice structure can be seen in 
Fig. ®. 



FIG. 4. Incommensurate phase. For k < n c the collinear 
phase (P) is established, and for larger values the continu- 
ous change of the wave vector, describing the intermediate 
phase, is plotted. The three-sublattice phase (J) is reached 
only asymptotically. The same route holds for the other two 
domains. 

Near point J the wave vector can be approximated by 

q c ~q + 1.32k(1, 0,0). (12) 

For systems with a dipolar strength as in CSN1F3, the 
transition to the incommensurate phase yields 



J' = 59 mK 



1. 



(13) 
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IV. EXCITATIONS 



B. Collinear antiferromagnetic phase 



Now we investigate the spin wave frequencies for the 
two commensurate phases within linear spin wave the- 
ory. The incommensurate phase resists such an analysis 
because of the infinite primitive cell. 

In linear spin wave theory the frequencies of the 
magnons can be derived via the the Holstein-Primakoff 
transformation (HP) which transforms the spin operators 
to Bose operators. Therefore, Bose operators a; and a]c3 
are introduced, which are given up to bilinear order by 

Sf = (S-a]ai), (14a) 

Sr = \[^(ai+4), (14b) 

Sf = -i^( ai -4) (14c) 

with the local spin vector S;. However, first the classical 
ground state must be determined, and then the the above 
expressions have to be inserted so that Sf corresponds to 
the local z-component of each spin. 



In the ground state there exist three domains (Fig. 
(|5|)). Spin wave theory is applied to domain A, where 
the spins arc orientated along the x-axis. The antiferro- 
magnetic modulation is given by the wave vector qi = 
-^(0, 1,0). For domains C and B the analogous wave 

vectors are q 2 = £ (1,1/V§,0) and q 3 = £ (1,-1/^3,0) 
respectively. 
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A. Ferromagnetic phase 

In this phase all spins are aligned ferromagnetically 
in the hexagonal plane. We choose the ^-component of 
the spins along the x-axis. After HP transformation the 
Hamilton operator has the form 

H = E fm + A * a l a * + \ B * (°q a -q + a \ a -J ( 15 ) 



with the coefficients 

A q = SA + 2S(J - J q ) + 2S(4 - J q ) + 

S{2A XX - A y y - A\ z ) (16) 
Bq = -SA + S(A% -A y y + 2iA^) (17) 



and the classical ground state energy 

Ef m = -NS 2 (J a + 4 + A™) 



(18) 



The planar anisotropy A does not contribute to this en- 
ergy. The spin wave frequency following Eq. (|l^), given 
by 



Eq = JAl - |B q p . 



(19) 



Due to the rotation symmetry in the plane the magnon 
frequency vanishes at the zone center. This spectrum 
becomes unstable for increasing exchange energy J' at 
point P which corresponds to a transition to the collinear 
antiferromagnetic phase (s. section III). 



FIG. 5. The ground state for CSN1F3 in the hexagonal 
plane is one of the three (A-C) shown configuration called 
domains. In domain A the two primitive vectors are repre- 
sented. The antiferromagnetic modulation can be described 
by qi, q2 and q3 for domain A, B and C respectively. 

The Holstein-Primakoff transformation can be applied 
for the crystallographic Brillouin zoneE3 by using the fac- 
tor e lqiX ' for the antiferromagnetic modulation 
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with the coefficients 
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(20) 



Aq — S(2Jq 1 — Jq J, 



q+qiy 



SA 



s(2 j qi - j q - j q+qi ) + s(2A% - a™ - A% + qi ) (21) 

Bq = S(Jq + q 1 — Jq) — SA + 

■^(^q+qi — Jq) + $ '(^q+qj — ^q) 

(22) 
(23) 



Cq — iSAq 

Dq = %SA» 



and the ground state energy 

^5 2 (j + j; 1 + ^). 



E* 



(24) 



Diagonalization of the quadratic Hamiltonian (Eq. 
(pp|)) can be_done by a generalized Bogoljubov 
transformatiorci The magnon frequency is given by 
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= ^i(n 1± n 2 ) 



(25) 



V. COLLINEAR ANTIFERROMAGNETIC 
PHASE WITH H / 



with 

fil = A q - S q + ^l q+qi - -Sq+q! + 8C q C q+qi 

and 

^ - (K - B$ - A q+qi + i^ +qi ) 2 + 

16[C q+qi (yl q+qi - -Bq+qJ - C q (Aq - Bq)] 
X [Cq (A q+qi + Bq+qJ - C q+qi (A q + Bq)} . 

Here we changed to the magnetic Brillouin zone, which 
is half the crystallographic, and therefore two branches 
appear. Note that the dipole energy lifts the degener- 
acy for q z because of the off-diagonal dipole com- 
ponent A^ z . The splitting of the two branches is of the 
order of the dipole energy. For wave vectors in the plane 
the off-diagonal component vanishes and the two magnon 
branches are degenerate. For wave vectors in the hexago- 
nal plane (q z — 0) stability of the collinear phase requires 

A q >0, A q >|B q |, q= (q x ,q y ,0). 

From these inequalities we recover the upper bound for 
the exchange energy J' (Eq. (|l^)). Note that for J and A 
much larger than J' and the dipole energy the inequalities 
are independent 9$-iJ and A. In Fig. (@) the dispersion 
relation is plottedtHl for CsNiF3 in two different directions 
of the Brillouin zone. 



0.20 




FIG. 6. Dispersion relation for CsNiF3 (vanishing field) for 
two directions. There is a discontinuity at the zone center. 

Note the discontinuity at the zone center, which is a 
result of the nonanalytic behavior of the dipole-dipole 
interaction (s. appendix). 



In this section we study the ground state and the dis- 
persion relation for the collinear antiferromagnetic phase 
in a homogeneous external magnetic field. First we in- 
vestigate the Neel phase and find that it is stable in the 
presence of a longitudinal field (parallel to the spin ori- 
entation) up to a critical value Bq. A transverse field 
always leads to a reorientation of the spins, i.e. a canted 
structure. In a second section we study the paramagnetic 
phase; especially we derive the instability boundary for 
arbitrary field direction at which the system undergoes 
a transition to a canted structure. Thus we obtain the 
corresponding wave vectors characterizing the interme- 
diate phases. In the last section we derive the stability 
equations for a general two-sublattice structure and the 
conditions under which conventional spin-flop phases can 
occur. 



A. Neel phase 

We consider domain A with the magnetic field along 
the x-axis. The classical ground state energy (Eq. (p^)) 
is not affected by the magnetic field because of the alter- 
nation of the spins. In the bilinear part of the Hamilto- 
nian (Eq. (|20|)) only the coefficient Z3 q changes to 

D q = iSAl z + ig L ^#o • (26) 

This affects the two branches of the dispersion relation 
in the following way: 

Eg™ = ^(fii±fi 2 ) (27) 

with 

^i = ^q-B q +A q+qi -B q+qi +8Cq C q+qi +2(g L fi B H ) 2 
and 

n 2 2 = (Al-Bl-Al + ^+Bl +cli ) 2 + 

16[C q+qi (A q+qi - B q+qi ) - C q (Aq - Bq)} 
X [Cq (A q+qi + B q+qi ) - C q+qi (Aq + Bq)} 
+ A(g Lf l B H ) 2 ((Aq + Aq+q,) 2 - (Bq - B q+qi ) 2 ) . 

The effect of the magnetic field is to lower (lift) the en- 
ergy of the lower (upper) branch. The collinear ground 

(2) 

state is stable as long as the second mode Eq ' is posi- 
tive. The energy of this mode vanishes with increasing 
field strength at first at the boundary of the magnetic 
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Brillouin zone at q4 = —(1,0,0). This implies a critical 
field of 

g L ^ B H c Q = 25 v /(^-^» +qi )( J 4--^») . (28) 

Surprisingly, the critical field depends only on the dipolar 
energy because the antiferromagnetic exchange energy 
obeys J q4 = J q4+qi = JL ■ The wave vector q4 describes 
a four-sublattice system, i.e. the primitive cell of the 
magnetic lattice consists of four spins. For CSN1F3 Eq. 
( p8| ) yields for the critical field 

= 47mT . (29) 

Above this value the collinear structure becomes unstable 
and an intermediate state occurs, which is investigated 
in the next section. The critical value is proportional to 
the energy gap at wave vector q4 present without field 
which then is increasingly lowered by applying a field. 
This phenomenon is already known for systems with an 
easy-.axis anisotropy and a magnetic field parallel to the 
spinsEa. 

B. Paramagnetic Phase 

Before studying the noncollinear phase, i.e. the spin- 
flop-like phase, it is instructive to investigate the param- 
agnetic phase. This phase is established for strong mag- 
netic fields in the hexagonal plane so that all spins align 
parallel to it. When the magnetic field strength is low- 
ered, the paramagnetic phase becomes unstable and a 
transition to a canted structure will occcur. From this 
instability point we obtain the wave vector characterizing 
the canted structure. Because the ground state of CSN1F3 
is not invariant with respect to a rotation around the 
spin-chain axis (recall that there are three domains A- 
C), the direction of the field plays a crucial role. Thus, we 
get a non-circular instability line, i.e. the absolute value 
of the critical field at which the paramagnetic phase be- 
comes unstable depends on the angle the magnetic field 
encloses with the x-axis. In consequence to the hexagonal 
symmetry it is sufficient to study a range of 60°. In the 
following tp denotes the angle between the external mag- 
netic field Ho and the x-axis. After a HP transformation 
the Hamiltonian reads: 

H = E PM + ^1 + 2 ( jB 1 a q + B q a q a -q) 

q 

(30) 

with the coefficients 

A q = 5(2( J - J q ) + A + 2 J - 2 J q ) + g L (M B H 
+S(2A XX - - sin 2 <pA™ - cos 2 ipA™ 
-sin2^A q ^) (31) 

B q = S(-A + A^ - sin 2 tpA* x - cos 2 tpA™ + sin 2tpA%> 
+2icos<M q z -2isin(M q 2 ), (32) 



and the classical ground state energy of the paramagnetic 
state 

Epm = ~NS 2 (J Q + J + A% x ) - g Lf i B NSH , (33) 

which is independent of the direction of the field. The 
dispersion relation for the paramagnetic phase is given 
by 

£ q = ^4|-|£ql 2 - (34) 

In Fig. (0) the excitation frequencies are shown for dif- 
ferent field directions (tp = 0). 




FIG. 7. Dispersion relation for CSMF3 in the paramagnetic 
phase. The field is aligned along the x-axis (tp = 0). The up- 
per curves are calculated for Ho = 2.1Hg x , the middle curves 
for Ho = l.Q6H§ x and the lower curves are just above the crit- 
ical field Hq*. Point K equals q = ^ (0, 0, 1). In graph (a) one 
can see the softening of the dispersion at an incommensurate 
wave vector. 

When the field is lowered the paramagnetic phase be- 
comes unstable below a critical field, given by the zero of 
the excitation energy (q z = 0, A q z = 0): 

(Jl^b Hq = 2S( J q - J - A xx + sin 2 <pA™ + cos 2 tpA™ 
-sin2tpA x v) . (35) 

Considered as a function of q the maximum of this ex- 
pression gives the critical value (dependent on the an- 
gle), and the wave vector q reflects the ordering of the 
phase below. Assuming a conventional spin-flop phase 
or an intermediate phase (two sublattice structure) we 
would expect the maximum value for the critical field at 
wave vectors describing the antiferromagnetic domains, 
i.e. qx, q2 or q3. However, the detailed analysis showed 
that the paramagnetic phase gets unstable at an incom- 
mensurate wave vector (\(tp) for certain field directions. 
Before considering the general case we study the special 
angles tp = 0° and ip = 90°. 



G 



tp = 0°: For this field direction Eq. (J35|) simplifies 



to: 



9Lf i B HS x = 2S(j' -4- a? + Ayy) . 



(36) 



This expression is proportional to defined in Eq. (Q). 
The maximum of this component is not achieved for qi 
but around qo. In the limit of vanishing dipole energy 
the incommensurate wave vector is given by (compare 
with Eq. ©) 



q(0°)~q - 1.32/s(l, 0,0), 



(37) 



i.e. the phase below the paramagnetic phase is an in- 
commensurate phase resembling a three sublattice struc- 
ture. Note that this incommensurate wave vector is dif- 
ferent from the wave vector found for the case of van- 
ishing dipole energy (note the minus sig in Eq. ( {L2[) ). 
For increasing dipole energy the wave vector (Eq. p7|)) 
changes continuously to q = 0, implying that for n —>■ 
the system orders ferromagnetically (s. III). 

For CsNiF3 the critical magnetic field can be evaluated 

to 



c. tp = arbitrary: Finally we turn to arbitrary an- 
gles, for which the situation turns out to be nontrivial. 
The complete dependence of the field direction q(y) for 
CsNiF 3 is plotted in Fig. §. 
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= 290mT 



at 



q(0) = -(1.023,0,0). (38) 
a 



The softening of this mode can be seen in Fig. (0(a)). 
This incommensurate wave vector happens to be near 
the wave vector q4, the four-sublattice structure, which 
describes an antiferromagnetic modulation along the x- 
axis. Note that q4 is the corresponding wave vector at 
which the Neel state becomes unstable (compare section 
A). 

b. tp = 90°: For magnetic fields parallel to the y-axis 
the critical value is given by 



9l^bH^ = 2S(J' -4- A™ 



A XX 

q 



(39) 



This expression is found to have its maximum value at qi , 
the antiferromagnetic wave vector of the collinear phase 
of domain A. For CsNiF 3 this yields 



Hq V = 340mT 



at 



q(90°) = qi 



(40) 



At Hq V the system undergoes a transition into a com- 
mensurate phase, precisely the two sublattice phase cor- 
responding to qi. When the antiferromagnetic exchange 
energy J' is lowered, the critical field vanishes for n < 
0.20. This limiting value agrees with the value for the 
transition between the ferromagnetic and the collinear 
antiferromagnetic ground state (Eq. (10)). In contrast 
to the previous case (tp = 0°), for this field direction there 
does not exist a region for n where an intermediate phase 
between ferromagnetic and paramagnetic phase occurs. 
Owing to the hexagonal symmetry the critical field for 
tp = 30° is the same as for tp = 90° (tp = —30°) but at 
the antiferromagnetic wave vector q2 (qa) characterizing 
domain C (B). 



FIG. 8. Wave vector q y and q x in the paramagnetic phase 
at which the phase gets unstable depending on the direction 
of the external magnetic field. We used the parameters for 
CsNiF3. tp measures the angle between the magnetic field 
and the a>axis. For tp > 7.8° the instability appears at point 
X. 

The wave vector components q x and q y are plotted as 
a function of the angle tp for the region in question. A 
critical angle tp c w 7.8° is observed above which the wave 
vector is locked to q2- Varying the angle tp between 0° 
and 7.8° the wave vector where the instability sets in 
changes continuously from q(0) to q(tp c ) = In Fig. 
@ the critical value for the magnetic field is plotted as 
a function of the angle for the same angular domain. 
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FIG. 9. Critical field below which the paramagnetic phase 
gets unstable for field direction between 0° and 30°. Note 
that there is no discontinuity at <p c = 7.8° but a slight kink. 
For smaller (larger) angles the paramagnetic phase changes 
to an incommensurate (commensurate) phase. The dashed 
curve indicates the critical field when assuming that the soft 
mode occurs at q2. 

Note that the critical field is a continuous curve even 
at the critical angle ip c . It has only a small kink at this 
point. The dashed curve results from Eq. ([35]) when 
it is assumed that the instability point occurs at q2 for 
the whole angular segment. Thus we can see that for 
ip < ip c the incommensurate structure is favored. Due to 
the inversion symmetry of the lattice together with q(<£>) 
there is a second modulation wave vector — q((f). The 
result of this investigation for a complete rotation of the 
magnetic field is summarized in Fig. ([l0|): The instabil- 
ity line of the paramagnetic phase is shown. There are 
angular regions where the paramagnetic structure under- 
goes a transition to a two-sublattice structure (annotated 
by the corresponding wave vector) separated by regions 
with incommensurate structures which are plotted with 
a thick line. The three dashed lines around zero indicate 
the Neel phases for a longitudinal field. 




q, 



p 





"Ax 




A/ x ' 



J 



FIG. 11. Angular dependence of the wave vector q(</?) on 
the instability curve from Fig. (10) for CsNiF3. Shown is the 
Brillouin zone from Fig. (6). The wave vector happens to be 
near the edge of the magnetic Brillouin zone for the special 
value of k, = 0.79. 

Now we want to derive the critical angle (p c , at which 
the commensurate and the incommensurate phases co- 
exist, for different strengths of the antifcrromagnetic ex- 
change energy J'. This value can be obtained by evalu- 
ation of the maximum of Eq. ( |35| ) near qi . Expanding 
the exchange energy and the dipole energy in a Taylor 
series we get 

4+qi " J ^+l J '(€-H) ( 41 ) 

^+ qi * A™ - G'(0.8W<f x + 0.112$) (42) 
4+ qi * A Z + G'(0.816f x - 0.112$) (43) 
* G'1.361^ , (44) 

where we introduced the dipolar strength G' = G/a^ 
(s. appendix). Inserting in Eq. ( |35| ) and setting the 
determinant of the Hess matrix zero we find the critical 
angle as a function of the exchange energy J' which is 
shown in Fig. (|l2|). 



FIG. 10. Angular dependence of the instability of the para- 
magnetic phase for CsNiF3-like systems. The thick-lined seg- 
ments on the instability curve correspond to the direction of 
the magnetic field Ho = Ho (cos </?, sin tp, 0) , at which the para- 
magnetic phase changes to an incommensurate phase. They 
are seperated by segments where a transition to a commensu- 
rate structure happens, especially a two-sublattice structure. 
The corresponding wave vector is given. The dashed lines 
denote the Neel state. 
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In Fig. ( |il| ) the full wave vector dependence within 
the Brillouin zone is shown for CSN1F3. Note that the 
wave vector is locked at the two-sublattice wave vectors, 
which cannot be seen in this representation. 
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FIG. 12. The dependence of the critical angle <p c on the 
angtiferromagnetic exchange energy J'. The region between 
the curve and <p = 90° denotes the commensurate structure 
right below the instability curve. 



of two sublattices on each of which the spin component 
along the field is the same (a = (3 in Fig. (13)). 



One can see that the angular segments of the commen- 
surate structure decrease when the antiferromagnetic ex- 
change energy takes values near the upper or the lower 
bound of the collinear antiferromagnetic phase. For 
CSN1F3 it happens that the commensurate segments have 
nearly the largest possible size. 



C. Spin-flop phase, intermediate phase 




In this section we want to examine the noncollinear 
phase, i.e. the region within the instability curve of Fig. 
@. In the following we investigate a magnetic field 
along the x-axis. Concerning domains B and C it is clear 
that for finite field strength the spins reorientate in order 
to gain energy from the interaction with the magnetic 
field, while domain A remaining in a Neel state cannot 
gain Zeeman energy. Thus the Neel state cannot be a 
stable state. In section A we showed via stability calcu- 
lations that a Neel structure is stabilized up to a critical 
value (Eq. (p8|)). Thus the Neel state for domain A 
cannot be a stable state but must be metastable. There 
is an activation energy necessary to flip domain A in a 
ground state emerging for domains B or C in the pres- 
ence of a field, due to the broken rotation symmetry in 
the hexagonal plane. The situation for a virgin probe is 
as follows: without magnetic field the system might be 
built up of the three domains with equal volume. Rais- 
ing the magnetic field (parallel to the spins in domain 
A) does not change the spin orientation in domain A but 
leads to a slight reorientation in domains B and C. Above 
the the critical value (Eq. (|2S|)) the spins in domain A 
flip (first order phase transition) to an orientation identi- 
cal to either domain B or domain C. For strong magnetic 
fields one finally enters the paramagnetic phase. When 
the magnetic field is decreased thereafter the spins or- 
der again in the two domains, but domain A is never 
formed again because of the metastability of this domain 
with respect to fields along the x-direction. The sys- 
tem ends up in a state where only domains B and C are 
present. In this scenario we neglected thermal effects, 
possible domain wall energies and defects in the crystal 
which should be taken into account for a proper descrip- 
tion. In fact, the experiments do not show this scenario; 
the system always ends in a three-domain state after a 
magnetic cycleo. 

For general field direction we expect from the investi- 
gation of the instability line of the paramagnetic phase 
that in the intermediate phase there are angular regions 
with incommensurate modulations separated by region 
with a commensurate two-sublattice modulation. Here 
we want to find out the conditions for existence of a con- 
ventional spin-flop phase, i.e. the ground state consists 



FIG. 13. Coordinate system for a general two-sublattice 
spin orientation. The spins of sublattice C\ (£2) 
enclose an angle of a (/3) with the magnetic field 
Ho = Ho (cos <p, simp, 0). 

Due to the large planar anisotropy A the spins reori- 
entate only within the plane, i.e. we need only one angle 
relative to the magnetic field for each spin. A two sub- 
lattice structure is characterized by the wave vectors qi , 
q 2 and q 3 . Introducing (Fig. (|3|)) the angle enclosing 
the x-axis and the magnetic field Ho by cp and the angle 
enclosing the field and the orientation of the spin on the 
first (second) sublattice L\ (£2) by a (/?), the ground 
state energy can be written in the form: 

Ef N = -NS 2 {J + cos 2 7 [J + A xx ] + 

sin 2 j[J^ + sin 2 {if + 8)A* X + cos 2 (ip + 5) A™ 

+ sin2(<p + 8)A^]} - g Lf i B NSH cos 7 cos 8 , (45) 

where q denotes one of the above mentioned two- 
sublattice wave vectors. Here we used the more appro- 
priate variables 7 and S, which are related to the original 
angles by 

a + (3 a — (3 



The minimum value of the ground state energy requires 
(d 7 Ef N = and d 5 E% = G): 

g L fi B H cos 8 = 2S[J^ - J' Q - A xx + sin 2 (tp + 8)A™ 

+ cos 2 ft? + 8)A™ + sin2((^ + S)A x i y ] cos 7 (46) 
gLVEiH cos 7 sin 8 = 5" sin 2 j[sm.2(ip + 8)(A x l x - A™) 

+2cos2(^ + 5)A^]. (47) 

These equations will be discussed now for a spin-flop and 
an intermediate spin structure. 

1. Spin- flop phase 

For a conventional spin-flop phase (a = /3 = 7 and 
5 = 0) these equations can be simplified 







gL^sHo 



4 



+ sm2ipA xy 
= sin2 ¥ >(A;f 



A XX 



q jC0S7 



Al v ) 



sin 2 <pA x * + cos 2 tpA y y 
(48) 

-2cos2^A q a . (49) 



The second relation gives a condition for the existence 
of a conventional spin-flop phase. By using the values of 
the dipole energies and the dipole relations given in the 
appendix we find the following solutions for Eq. (M9): 





qi q2 q3 


stable 
unstable 


TT 67T TT (TT bTT JLJ.7T 
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In the second row the angles are given for the correspond- 
ing wave vector. These solutions represent magnetic field 
directions perpendicular to the spin orientation for van- 
ishing field. In the third row we tabulated the solutions 
which correspond to field directions along the spin orien- 
tation. The latter solutions represent no minimal ground 
state energy for exchange energies J' < 26m K. This can 
be seen from Eq. (^8|), where the term in brackets then 
is negative. Even for Hq = and e.g. cp = 0° the ground 
state energy 



E 9 



-NS 2 {J + J^+A™) 



is larger than the ground state energy for the Neel state 
(Eq. (0)). Thus, there is no spin-flop phase for longitu- 
dinal fields in CsNiF 3 . 

Inserting the values (stable) for the angles from the 
table in Eq. ( p§|) we get the static equilibrium relation 
between the magnetic field and the angle 7 



9lVbHq 



2S[J^-4-A* * + A™} cos 7 . 



(50) 



Note that this relation holds for all stability values of 
ip due to the hexagonal symmetry, as can be seen by 
applying the dipole relations from the appendix. So we 
conclude that for most angles ip there does not exist a 
spin-flop phase except for the above found solutions. 



2. Intermediate phase 

In the general case of two independent angles we have 
to solve equations p§) and (^) . As a result one obtains 
the dependence of the angles a(Ho, ip) and ft(Ho, ip) upon 
the field strength and the direction. Since the solution 
for the general case is nontrivial we want to study only 
two special situations, namely ip — 0, ip = 90° and a 
modulation wave vector qi for domain A. In these two 
cases the equations of stability (Eqs. (48) and (|49|)) can 
be solved analytically. For <p = 90° we find that the 



intermediate phase has a higher ground state energy than 
the conventional spin- flop phase (a = 0). For ip = the 
ground state energy of the intermediate phase exceeds the 
energy for the Neel state for vanishing fields, signifying 
the absence of such a state. The complete investigation 
of the general case is left for future studies. 



VI. SUMMARY 

We have studied the ground state for a planar, quasi 
one-dimensional ferromagnet with weak antiferromag- 
netic exchange and dipole-dipole interaction. For vanish- 
ing exchange energy the spins order ferromagnetically in 
the hexagonal plane. With increasing exchange energy 
the ground state passes through a collinear phase with 
three possible domains and switches over to an incom- 
mensurate phase, where the characterizing wave vector 
moves continuously from the points P, X or X' to the 
three-sublattice wave vector at point J (Fig. (||)). The 
three-domain structure of the collinear phase persists in 
the incommensurate phase. 

Via linear spin wave theory we obtain the dispersion 
relation for the ferromagnetic and the collinear antifcr- 
romagnetic phase (Hq = 0). The range for the collinear 
antifcrromagnetic phase for CsNiFalike systems is cal- 
culated from ground state energy calculations and from 
stability conditions concerning to excitations 



7mK < J' < 59mK , 



(51) 



which is consistent with J' = 25 mK found in recent 
measurements^ 2 !. The upper and lower bounds are inde- 
pendent of the ferromagnetic exchange energy J and the 
planar anisotropy energy A, as long as they are much 
larger than the dipole energy and the exchange energy 
J'. 

In the following we studied the magnetic phase dia- 
gram of the collinear antifcrromagnetic phase. For ho- 
mogeneous external magnetic fields parallel to the spin 
orientation the Neel phase is stable up to a critical field 
which depends only on the dipolar energy (Eq. (p8|)). 
The instability is shown by the softening of a mode with 
wave vector q4 which describes a transition to a four- 
sublattice structure. For magnetic fields transverse to 
the spins a reorientation sets in for all finite field values. 

By stability investigations of the paramagnetic phase 
we obtained a non-circular instability line for fields in 
the hexagonal plane. There is a strong dependence on 
the field direction on the critical field and the structure 
of the phase below the paramagnetic phase. For certain 
angular domains the system changes in incommensurate 
structures which are separated by commensurate (two- 
sublattice) structures (Fig. (|Io|)). We showed that the 
commensurate angular domains are intermediate phases 
while conventional spin-flop phases exist only for mag- 
netic field directions <p = 30° ±-jn60°. This is in contrast 
to the result of Yamazaki et al£3, who found conventional 
spin-flop phases for all field directions, even for a longitu- 
dinal field where an incommensurate phase should set in 
from our calculations. This may result from their semi- 
classical model which does not consider the full nature of 
the dipole-dipole interaction. 

As a result, the magnetic phase diagram reveals an in- 
teresting structure of different phases. The investigation 
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of the complete magnetic phase diagram (intermediate 
phase) is left for future studies. 
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-£|!^( q + G r(q + G)% (■ 

G V 



2/3 
47T 



|q + G|' 



Here I refers to the direct lattice and G to the reciprocal 
lattice. The prime excludes the value I — and we used 
the Misra functions 



ip n (x) = I dtt n e~ xt 



Expanding the dipole tensor for small wave vectors 
yields: 



VII. DIPOLE-DIPOLE INTERACTION 

In this section we summarize the most important re- 
lations concerning the dipole-dipole interaction needed 
in the analysis. We consider a simple hexagonal lattice, 
with lattice constant a in the triangular plane and lattice 
constant c perpendicular to the plane. The direct lattice 
is paramatrized by 

x, = o(li,0,0) + |(Zi,V^2,0) + c(0,0,I 3 ), 

h = 0,±1,±2,... 
and the reciprocal lattice by 

2n 2n 2ir 

G = — (hi, 0,0) + ^r(2h 2 - h u 0) + —(0, 0, h 3 ), 
a V3a c 

hi = 0,±1,±2,... 
The Fourier transform of the dipole energy (Eq. (0)) 

<] 



G 



with 



lira J ( V e iqx ' ( — 

-"Hi V |x 



5 af} 3afxf 



X;| 3 Ix-X/I 5 



(52) 



G = (gh^B) 2 

is found by the Ewald summation techniqueEl. There 
the summation is split into a part over the direct lattice 
and a part over the reciprocal lattice and, by usins-in an 
intermediate step the generalized theta functiorJ13, we 
derive the following relation (03 = -^a 2 c): 

a' , XaJ J 



A « G 



-2tt 



27F " a I 77 I I 



7T 



a Za0 \ ' 



2/ 
i 3 I 



J2G a G P Vo f^j^) + <%V) (53) 



-2tt 



q a q §a 
~ 2 3~ 



aiat ■ 



(54) 



where we have introduced a lattice dependent term ai a t- 
This expression is non-analytic due to the semiconver- 
gence of the dipole sum in three dimensions. The value 
at q = depends on the direction of the limiting process 
expressing the shape dependence of the system. For a cu- 
bic lattice the lattice dependent term vanishes a; at = 0, 
so the dipole energy for a ferromagnetic structure is zero. 
In contrast, this term does not vanish for a hexagonal lat- 
tice. For non-spherical surfaces this equation has to be 
corrected by the demagnetization factor. 

In the following table the dipole components are listed 
for selected wave vectors in units of G/a 3 and for k = 
a/c = 2.39, the ratio of the lattice constants for CsNiF 
(«. = 0) 





q = 


qo 


qi 


q2 


q3 


q4 


Axx 


-2.805 


-5.946 


-3.392 


-7.223 


-7.223 


-7.520 


A vv 
q 


-2.805 


-5.946 


-8.501 


-4.669 


-4.669 


-4.373 




9.798 


11.893 


11.893 


11.893 


11.893 


11.893 


A*y 











-2.212 


2.212 






Due to the sixfold rotation symmetry of the lattice, we 
can derive a general relation between dipole components 
in different directions: 



A AVV _ 0411 

q(v) ~ dyi q(v+60°) 



q(v+60°) 



(55) 



q(y) denotes a wave vector along direction tp in the Bril- 
louin zone. From this formula we obtain relations for 
different points in the Brillouin zone: 



11 



4A» = 3^ + A™ - 2y/3A% (56) 

4A xx = + 3j4 to _ 2 V3A^ (57) 

AA xx = A xx + 3A yy + 2v ^» ( 58 ) 

= 3(A^ - A^) + 2V3A% . (59) 
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